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ABSTRACT 

Purely photometric models can be used to determine the binary parameters of eclips- 
ing cataclysmic variables with a high degree of precision. However, the photometric 
method relies on a number of assumptions, and to date there have been very few 
independent checks of this method in the literature. We present time-resolved spec- 
troscopy of the P = 90.9 min eclipsing cataclysmic variable OY Carinae obtained 
with X-shooter on the VLT, in which we detect the donor star from K i lines in the 
J-band. We measure the radial velocity amplitude of the donor star K2 — 470.0 ± 2.7 
km/s, consistent with predictions based upon the photometric method (470±7km/s). 
Additionally, the spectra obtained in the UVB arm of X-shooter show a series of Fe 1 
and Fe 11 lines with a phase and velocity consistent with an origin in the accretion 
disc. This is the first unambiguous detection at optical wavelengths of the 'iron cur- 
tain' of disc material which has been previously reported to veil the white dwarf in 
this system. The velocities of these lines do not track the white dwarf, reflecting a 
distortion of the outer disc that we see also in Doppler images. This is evidence for 
considerable radial motion in the outer disk, at up to 90km/s towards and away from 
the white dwarf. 
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1 INTRODUCTION 

The study of close binary stars is driven by the need to 
understand the numbers and properties of the many pos- 
sible outcomes of binary star evolution. This is of interest 
in many astrophysical contexts, such as the need to deter- 
mine the properties of Type la supernovae progenitors as 
a function of redshift, and the detection of gravitational 
waves, for which clos e binaries are predicte d to be among 
the strongest sources l|Nelemans et alj|2004l ) . A useful pop - 
ulation are the cataclysmic variables (CVs: I Warner! \l 995). 
These close binaries consist of a white dwarf accreting mat- 
ter from a main sequence donor star, and provide us with 
a large, homogeneous and easily observed population which 
can be compared to theoretical predictions. 

The binary parameters of eclipsing systems can be 
determin ed with high precision vi a a purely photometric 
method dBailevI I l979l ; ISmakl 1 19791 : ICook fc Warnerl 1 1984 
I Wood et al.lll98E ) in which the eclipse features in the light 
curve are used to infer the geometry of the system. The num- 
ber of known eclipsing systems has increased significantly in 



recent years thanks to programmes such a s the Sloan Dig- 
ital Sky Survey (SPSS: ISzkodv et alj|2009h . and so greater 
relia nce has been placed o n this photometric technique (see , 
e.g.. iLittlefair et alj |2008| ; ICopperwheat et all |2010| . 1201 ll ; 



ISavourv et alj 2011 ) . The photometric method is based on 
two key assumptions: first, it is assumed that we see the bare 
white dwarf and can therefore measure its radius, and hence 
its mass, using models of white dwarf mass-radius relations. 
Secondly, it is assumed that the bright-spot where the gas 
stream hits the accretion disc lies on the ballistic trajectory 
of the gas stream, which follows a path set by the binary 
mass ratio. To date there have bee n very few independen t 
checks of the photometric method. IWade fc Hornd (Il988l ) 
presented a radial velocity study of the secondary star in 
Z Cha, and found a value for the donor star orbital velocity 
K2 = 430±16km/s . In comparison, the photometric study of 
IWood et all l|l986l ) found K 2 = 371.0 ± 2.7km/s, a discrep- 
ancy of more than 3a. This was partly a ttributed to an over - 
estimation of the white dwarf radius in IWood et all (|l986t) . 
but later work accounting for this (increasing the photomet- 
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ric determination of K2 to 389 - 406k m/s) was still ^2cr di f- 
ferent from the spectroscopic finding l|Wood fc Hornc 1990). 
More recent spectroscopi c studies (IP Peg, Watson et al.l 
120031 ; CTCV J1300-3052, ISavourv et afll201lh have shown 
good agreement with photometric results. However, the 
small number of test cases to date is concerning, given the 
ubiquity of the photometric method. There have been some 
attempts to verify the method by inferring the white dwarf 
orbital v elocity (Ki) from th e velocities of the emission lines 
(see, e.g. lTulloch et ai1l2009h . However, these measurements 
are generally based on the assumption that the accretion 
disc is perfectly symmetric around the white dwarf. This is 
often not the case, and while there have been a ttempts to 
correct for a non-symmetrical disc (|Marsh| [l988) we would 
still consider these K\ determinations to be less reliable than 
measurements of the donor star. 

The capabilities of modern spectrographs far exceed 
the instrument used by IWade fc Hornd (|l988l ). In partic- 
ular, the wavelength range and throughput of X-shooter 
mounted on the ESO Very Large Telescope (VLT) enables 
us to determine a value of K2 with a precision comparable 
to that which can be currently achieved via the photometric 
method. We therefore chose to measure this orbital velocity 
for the donor star in a CV and thus provide a more precise 
test of the light curve fitting method. The system we selected 
was OY Carinae (henceforth OY Car), a bright and short- 
period (P = 90.9 min) eclipsing system with a prominent 
white dwarf eclipse. The best photometric determinations 
fo r the binary parameter s in OY Car to date were published 
in lLittlefair et all (|2008l ). 



2 OBSERVATIONS AND REDUCTION 

We observed OY Car for a total of 1.8h on the night o f 
2010 February 10 with X-shooter l|D'Odorico et al .11 20061 ) 
mounted on VLT UT2 (Kueyen). The X-shooter spectro- 
graph consists of 3 independent arms (UVB, VIS and NIR), 
giving a simultaneous wavelength coverage from 3000 to 
24800A. We obtained 28 consecutive spectra of 185, 180 and 
205s in length for the UVB, VIS and NIR arms respectively. 
The cycle time per exposure was ~235s. We binned by a 
factor of two both spatially and in the dispersion direction 
for the UVB and VIS arms, and used slit widths of 1.0", 
1.2" and 0.9" for the UVB, VIS and NIR arms respectively, 
resulting in FWHM spectral resolutions of 0.7, 1.0 and 2.9A. 
Weather conditions were excellent, with sub-arcsecond see- 
ing and photometric transparency throughout. 

We reduced these data using version 1.3.7 of the X- 
shooter pipeline. The standard recipes were used to opti- 
mally extract and wavelength calibrate each spectrum. We 
removed the instrumental response using observations of the 
spectrophotometric standard star CD 71. One complication 
was that, owing to an oversight, the data were obtained in 
'stare' mode, rather than by nodding the object along the 
slit, as is typical for long slit infrared spectroscopy. The re- 
sult is the sky subtraction is poorer than is usually possible 
with X-shooter. We divided our spectra by the spectrum of 
a telluric standard star taken at a similar airmass to reduce 
telluric absorption. Fortunately, most of the spectral lines we 
use are displaced from the worst telluric features; the optical 



data are not affected. We plot the averaged and continuum 
subtracted spectra from all three arms in Figure [T] 



3 RESULTS 

3.1 The orbital velocity of the donor star 

In order to measure K 2 in Z Cha, IWade fc Hornel (|l988l ) 
used the AA 8183, 8194 A Na I doublet to track the motion 
of the donor star. The wavelength coverage and sensitiv- 
ity of X-shooter means other absorption lines are available, 
and infrared spectra of M, L and T dwarfs show promising 
features in the J- band (|Cushing et al.l [20051 . table 6) which 
have been used in previou s CV studies to detect the donor 
star (|Littlefair et al.ll2000l ). We find the K I doublets at AA 
11690, 11772 A and AA 12435, 12522 A to be much stronger 
features in our OY Car data compared to the Na I doublet, 
and we used these features to track the donor star. We plot 
the trailed spectra around these absorption lines in Figure 
[2] Some residual telluric contamination is apparent in these 
trailed spectra: the 12435A line is particularly affected. We 
masked out these sections of data when applying our model 
fits. 

We fitted the spectra with a model consisting of a 
straight line and Gaussians for the K absorption lines, al- 
lowing the position of each Gaussian to change velocity ac- 
cording to 

V = 7 + K 2 sm(2ir(p), 

where 7 is the offset of the line from its rest wavelength 
and <f> is the orbital phase of the spectrum. We obtained the 
rest wavelengths from the NIST Atomic Spectra Database 
ijRalchenko et all |2010| ). The spectra were fitted si multa- 
neou sly using a Levenberg-Marquardt minimisation (|Pressl 
2002), and the results are listed in Table[T] The zero point for 
the orbital phase was a free parameter in these fits, and we 
find it to be at HJD = 2443993.54892(4). We used the orbital 
eriod from the quadratic ephemeris given in lGreenhill et ail 
2006). We give the results for each line fitted individually 
with a single Gaussian model, as well as a combined fit in 
which the four lines are fitted together with a common 7, 
K2 and velocity width. The combined fit gives a value of 
K2 = 470.1 ± 2.0km/s for the orbital velocity of the donor 
star. The model fits typically give a reduced \ 2 of ~1.4, so 
we rescale the uncertainties to give a reduced x 2 of 1, and 
it is these values we quote in the table. We also make a sec- 
ond determination of the combi n ed fit, for which we use th e 
bootstrap method jEfronl Il979l ; lEfron fc Tibshiranil 1 19931 ). 
In this method we determine the parameter values and un- 
certainties via fitting a large number of individual sets of 
spectra which have been resampled from the original set. 
The two combined fits give parameter values which are in 
excellent agreement, with the uncertainties slightly higher 
for results determined via the bootstrap method. We used 
the combined fit as a starting point to fit each spectrum 
individually, and plot the resultant radial velocity curve in 
Figure [3] along with the sinusoid corresponding to our orbital 
solution. The sinusoid is a good fit to the radial velocity - 
there is little to no evide nce for irradiation as i s observed in 
other CVs (e.g. EX Dra. iBillington et al.lll996l ). 

The combined (bootstrap) fit gives a value of K2 = 
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Figure 1. Averaged and continuum subtracted spectra. The top, middle and bottom panels show the spectrum obtained in the UVB, 
VIS and NIR arms of the spectrograph, respectively. We give the wavelength in Angstroms, and the panels arc of different widths so 
as to keep the scale consistent. We identify the most prominent features, and mask out the regions of the spectrum which are most 
contaminated by telluric features. 



470.0 ± 2.7km/s for the orbital velocity of the donor star. 
When we fit the lines individually we find consistent values 
for K2 for three of the four lines, with the 11690A value 
around 1.5a lower than the value for the combined fit. Our 
determination of the systemic velocity 7 from the combined 
fit is 66.2 ± 2.5km/s, although this determination is less se- 
cure since three of individual fits result in 7 values which 
differ significantly from the combined value. The formal er- 
rors we quote would seem to be an underestimation of the 



true uncertainty in this parameter at least, probably due to 
the telluric absorption. 

I Wood fc Hornd jl990T ) showed that the photometrically 
determined value of K2 depends upon the white dwarf ra- 
dius, the determination of which depends on the limb dark- 
ening coefficient used. For OY Car they reported K2 values 
of 470 ± 8 and 455 ± 9km/s for a limb darkening coeffi- 
cient of uwd = or 1 respectively, where the white dwarf 
intensity distribution 1(8) = /i(1 + iiwd — uwd cos 8). How- 
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Figure 2. Trailed spectra obtained in the NIR arm. We plot the regions around the AA 11692, 11778A and AA 12437, 12528A K I 
doublets. These absorption lines move in phase with the M-dwarf donor. The scale is such that the continuum is white with darker 
absorption features. 'JV' denotes the number of the spectrum. The N = 19 spectrum was obtained around phase 0. 



Table 1. Parameter fits to the four K I lines observed in the 
NIR arm. We list the systemic velocity 7, the donor star orbital 
velocity K2, and the full- width half maximum (FWHM) of the 
absorption line, which we also express in velocity units. We list the 
results from fitting each line individually with a single Gaussian 
model. The uncertainties quoted are the formal errors, scaled to 
give a reduced x 2 of 1. We also list the results of combined fits 
obtained with a model consisting of four Gaussians. We give two 
combined results: one in which the uncertainties are the scaled 
formal errors, as for the individual fits, and a second in which we 
determined the parameter value and uncertainty via a bootstrap 
method. The bootstrap method gives slightly larger uncertainties. 



Wavelength 


7 


K 2 


FWHM 


(A) 


(km/s) 


(km/s) 


(km/s) 


11690.22 


74.5 ±4.6 


457.1 ±6.2 


235.6 ±10.8 


11772.84 


58.8 ±2.7 


473.2 ±3.4 


246.6 ±6.4 


12432.27 


77.4 ±4.3 


464.1 ±5.6 


243.7 ±8.5 


12522.14 


67.4 ±3.0 


470.5 ±4.1 


206.9 ±6.7 


Combined (x 2= l) 


66.2±1.6 


470.1 ±2.0 


238.9 ±3.8 


Combined (Bootstrap) 


66.2±2.5 


470.0 ±2.7 


239.1 ±4.4 



ever, they used a zero temperature mass-radius relationship 
for the white d warf, which introduced a significant bias. 
iLittlefair et al.l (|2008l ) revisit ed the system parame ters of 
OY Car, and improved on the lWood fc Hornel <jl990l ) deter- 
mination by using a more realistic finite temperature model 
for the whi te dwarf. They ass umed uwd = 0.5 and Two = 
16, 5000K (jHorne et al J 1 19941 ) and find Ki = 470 ± 7km/s. 
This is in excellent agreement with our spectroscopic value 
of 470.0 ± 2.7km/s. 



3.2 Rotational broadening of the donor star 

As an additional test of the photometric parameters, we 
measured the rotational broade ning of the d o nor s tar V ro t 
following the method detailed in lMarsh et alj (Il994t ). Using 
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Figure 3. Radial velocity curve for the secondary star. The solid 
line is the sinusoid corresponding to our combined fit orbital so- 
lution, as determined in Section |3. II 



the model fit detailed in Section [XT] we applied an offset to 
each spectrum to remove the orbital variation and then av- 
eraged the spectra. We then applied an artificial rotational 
broadening to the spectrum of a template star, and sub- 
tracted the result, multiplied by a constant representing the 
fraction of flux in the donor star, from the averaged OY 
Car spectrum. The constant was varied to optimise the sub- 
traction, and the \ 2 difference between the residual spec- 
trum and a smoothed version of itself is computed. This 
process was repeated for a range of rotational broadenings. 
By compari ng our spectra with t he spectroscopic sequence 
presented in lCushing et all |2005l ) we identify the donor in 
OY Car to be a late M-dwarf (M8/M9) or possibly an early 
T-dwarf. We rule out an earlier type M-dwarf since in these 
stars the K I features are weaker, and there are a num- 
ber of Fe absorption lines of comparable strength around 
the AA 11690, 11772 A doublet which we do not observe 
in our OY Car spectra. We therefore used used LHS 2021 
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Figure 4. Curves of reduced \ 2 versus V ro t sin i, obtained via 
the rotational broadening measurement procedure described in 
Section l3,2l For the template star we use LHS 2021, an M8 dwarf 
(solid curve) and LHS 2065, an M9 dwarf (dashed curve). We 
assume a limb darkening coefficient of 0.6. The s olid vertical line 
marks the value of v sin i determined using the IWood fe Horn j 
l ll990h binary parameters. The vertical dashed lines m ark the la 
error on this value, using the iLittlefair et al parameter 
uncertainties. The vertical dot-dashed lines show the uncertainty 
when we use our tighter constraint on Kn. 



and LHS 2065 as our template stars, which are of spectral 
type M8 and M9 respectively. An additional input p arame- 
ter is t he limb darkening coefficient for the donor star. IClaretl 
( 1998) listed limb darkening coefficients for low mass stars, 
and from these tables we find a value of it2 ~ 0.6 - 0.7 is 
an appropriate J-band coefficient for the donor. We plot in 
Figure [4] the resulting \ 2 curves for U2 = 0.6. The minima 
of these curves give our preferred values of V ro tSaxi, which 
we find to be HOkm/s and 109km/s for the M8 and M9 
templates, respectively. 

The quantity V ro t sin i is related to the mass ratio q and 
Ki via 

R 2 

V r ot sin i — K 2 [q + 1) — 



l|Horne et al.lll986l ), where for a Ro che-lobe filling star R2/a 
is a function of q (jEggletonl Il983j). Using the photom etri- 
cally determined parameters from Littlefair et alj (|2008l ). we 
find Vrot sin i = 108 ± 6km/s. This uncertainty decreases to 
±5km/s when we use our spectroscopically determined con- 
straint on K-2. We indicate this photometric determination 
of Vrot sin i in Figure [4] The figure demonstrates that, for 
both template stars, our spectroscopically determined value 
of Vrot sin i for U2 = 0.6 is consistent with the photometric 
value to within the uncertainties. We should note that this 
U2 is appropriate for the continuum of the donor star in the 
J-band, but the coefficient for the absorption lines may be 
somewhat different, which would change the optimum spec- 
troscopic value of Vrot sin i. However, we find changing 112 
in increments of 0.1 shifts the optimum value by only 1 - 
2km/s, and so given the uncertainty in the photometric de- 
termination of V r ot sin i, the spectroscopic value will still be 
consistent even if the coefficient for the lines is very different 
to our assumed value. 



3.3 Doppler maps: emission from the disc and 
donor star 

We cre ate Doppler maps using a Maximum Entropy Method 
(MEM. iMarsfi fc Home! Il9881 ; iMarshl l200ll ) for the 4686A 
He II emission line, and the Ca II triplet at AA 8498, 8542, 
8662A. We plot the trailed spectra and Doppler maps in Fig- 
ure [5] The Ca II map is a single, common map to represent 
all three lines, with a relative scaling of 1.00, 1.25 and 1.06 
for the three lines respectively. If we examine the Ca II map 
first, we see a ring of emission which originates in the accre- 
tion disc. The velocity of this emission is ~900km/s, which is 
towards the high end of the range observ ed in other CVs (e.g. 
600km/s in U Gem, Marsh et al.l 19901 ; see also many other 
examples in lKaitchuck et al.lll994l ). This implies a small ac- 
cretion disc, and is consistent with the short orbital period 
and large white dwarf mass of OY Car. We also see the 
'bright spot' of emission where the gas stream impacts the 
accretion disc, as well as Ca II emission from the donor star. 
This donor emission can be identified as a sharp S-wave com- 
ponent in the trailed spectra, and in the Doppler map we see 
that it is concentrated in the irradiated face of the star. We 
also created (but do not plot) Doppler maps of the Balmer 
lines. These show the same features which we observe in the 
Ca II map, but the donor emission is not as sharply defined. 

Turning to the He II emission, in the trailed spectra we 
see an S-wave component with an amplitude and phase that 
is inconsistent with a donor star origin. In the Doppler map 
this emission resolves to a single component in the bright 
spot. Both the Ca II and He II bright spot emission have 
kinematics which place them between the disc and stream 
kinematics. The emission in both maps is quite smeared out, 
but the Ca II emission may be slightly closer to the velocity 
of the incoming gas stream material. This would suggest 
that the emission in He II is from the mixture of gas stream 
and disc material further downstrea m, as is seen in other 
CVs (e.g. U Gem. lMarsh et al.lll990l ). 



3.4 The 'iron curtain' in the accretion disc 

In the UVB arm (3000 - 5600 A), we observe many absorp- 
tion lines with a common phase and velocity. The velocity 
is such that these features cannot be attributed to either 
the white dwarf or the donor star, and so we conclude they 
originate in the accretion disc. The majority of these fea- 
tures lie in a wavelength range of 3100 - 3500A (Figure©, 
but additional lines can be observed longwards of this to 
~5000A. 

In Table [2] we list our identifications for these absorp- 
tion lines. We find the majority to be Fe I and Fe II lines. 
These identifications, coupled with the evidence for an ac- 
cretion disc origin, suggests we a re seeing the ' iron c urtain' 
first identifed in this system by iHorne et alj i| 19941 ). They 
obtained ultraviolet observations of OY Car with HST, and 
found that the emission from the white dwarf was veiled by a 
forest of Fe II features, which they attributed to absorption 
by material from the outer accretion disc. This effect has 
never been directly detected before at optical wavelengths, 
and our optical detection is at a much higher spectral reso- 
lution than the ultraviolet HST data. At these longer wave- 
lengths we find the spectra t o be dominated by Fe I features, 
whereas IHorne et all |l994h found the majority of the ab- 
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Figure 5. Top: trailed spectra around the 4686A He II emission line (left) and the AA 8498, 8542, 8662ACa II triplet (right). The 4713A 
Hel line is also visible in the He II trail. The scale is reversed for these plots compared to the other trailed spectra in this paper, with 
the dark lines showing emission features rather than absorption. 'N' denotes the number of the spectrum. The N = 19 spectrum was 
obtained around phase 0. Bottom: Doppler maps for the He II (left) and Ca II (right) emission lines. In the Ca II we produce a combined 
map from the three lines of the triplet. The 'spoking' evident in this map is an artefact due to the relatively low phase resolution of 
our spectra. We overp lot in blue the Roche lobe of the donor star and the streams for a model with q = 0.102, K\ = 48 and A*2 = 470 
l|Wood fc H ome 1990). The streams denote the velocity of material in the gas stream itself (lower) and the velocity of the disc along the 
stream (upper). We mark with a cross the centre of mass, and with a square the position of the white dwarf. Points of equal distance 
along the two streams are marked, in units of the distance to the inner Lagrangian point from the centre of mass of the white dwarf. 
The small and large points mark increments of 0.02 and 0.1 of this distance, respectively. 



sorption features in their data to be due to Fe II lines. The 
same effect could explain the absorption cores in the Balmer 
lines, which are often observed to drop below the level of the 
continuum in eclipsing CVs. 

Using the rest wavelengths listed in Tabled we fit each 
spectrum individually in the same way as we did for the 
donor star in Section 13.11 We plot the resulting velocities 
in Figure [7) The sinusoid we plot for comparison in this 



figure has an amplitude of 90km/s and a systemic velocity 
consistent with that derived for the donor star, and is an 
adequate fit to most of the velocities. There is a significant 
departure from this sinusoidal behaviour for the few points 
prior to the eclipse. This is in line with a disc origin for these 
absorption features since at these phases the line-of-sight 
will be through the bright spot, and so our measurements 
will be distorted by the complex dynamics of the disc-stream 
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Figure 6. Bottom: Trailed spectra obtained in the UVB arm between 3100 and 3500A. We observe a series of absorption lines, with a 
phase and velocity consistent with an accretion disc origin. The scale is such that the continuum is white with darker absorption features. 
'N' denotes the number of the spectrum. The N = 19 spectrum was obtained around phase 0. Top: Averaged spectrum over the same 



interaction region. As well as the velocity, we note that the 
lines seem to vary in strength with phase: Figure [6] shows 
them to be fainter (and also slightly broader) in the few 
spectra leading up to phase 0.5 (N = 8) and phase (N = 
19), but much stronger in the intermediate phases. 

After removing the orbital velocity shifts and averag- 
ing the spectra, we fitted a subset of the iron lines and 
determined their full-width half maximum (FWHM) to be 
85.5 ± 2.9 km/s. For this subset we excluded lines which 
were potentially blended. We find no significant difference 
between the FWHM of the lines we identify as Fe I com- 
pared with those we identify as Fe II. 



3.5 A simple model for the iron curtain lines 

We computed a simple model which determines the absorp- 
tion for a thin gas layer which rotates with the disc. This gas 
layer is assumed to be uniform in height, and is split into 
multiple sub-layers so the effect of shear ca n be accounted 
for. W e use the binary parameters given in iLittlefair et all 
(2008) for these calculations. We considered various loca- 
tions for the absorbing gas layer. Firstly, we assumed it was 
very close to the white dwarf, and in the plane of the disc 
itself. In this configuration, the light which is absorbed orig- 
inates in the lower hemisphere of the white dwarf and passes 
through the region of the (inclined) accretion disc which is 
on the line-of-sight to the observer. However, in this config- 
uration our model predicts much broader absorption lines 
than we observe. In order for us to obtain line widths that 
match the observation, we require the observing region to 
be located in the outer accretion disc. Since the binary is in- 
clined, this also requires the absorbing region to be elevated 
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wavelength range, with the velocity variations subtracted. 
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Figure 7. Radial velocity curve for the absorption lines observed 
in the UVB arm. The solid red line is the best sinusoid fit. The 
dashed blue line is the r adial velocity curve fo r the white dwarf, 
assuming K\ = 48km/s l|Wood fe Hornelll990l l. 



8 CM. Copperwheat et al. 



Table 2. Absorption features identified in the spectra obtained in 
the UVB arm. The rest wa velengths are obtained from the NIST 
Atomic Spectra Database (Ralchcn ko et al1l201Clh . 
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Fe I 


3433.57 


Fe I 
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Fe I 


3459.91 
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3570.10 
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Fe I 


3758.23 


Fe I 


3761.41 
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4045.81 


Fe II 


4583.84 


Fe II 


4923.93 


Fe II 
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above the plane of the accretion disc, otherwise it would not 
intersect the light emitted from the white dwarf. In Figure 
[8] we plot a cartoon which illustrates the small region of the 
disc we are probing. We define a height h above the accre- 
tion disc for the gas layer. The h we select also determines 
a minimum distance rj„ for the observed region of the gas 
layer. Note that the gas layer is not necessarily truncated 
at ri n , but any material closer in does not lie on the line- 
of-sight to the white dwarf and so is not observed. We do 
truncate the gas layer at the outer radius of the accretion 
disc, which we assume to be ~20 white dwarf radii. Note 
also that as the binary phase changes the azimuth of the 
region we probe also changes, so if we obtain an entire orbit 
of time series spectra we map out an annulus of material 
around the disc. 

In Figure [9] we plot the region around three represen- 
tative iron absorption lines in the averaged spectrum. We 
choose lines which are strong and relatively isolated. Some 
of the lines show possible emission wings as well as the ab- 
sorption line: if this emission is real it presumably originates 
in the material which is not on the line-of-sight to the white 
dwarf. We also plot a series of model predictions. We exam- 
ine first the solid lines. For all of these models we assume an 
RMS line width of lOkm/s. We consider this an appropriate 
value to account for the turbulence in the disc. We plot the 
model in which the absorption is in the disc plane, as well as 



WHITE DWARF 




Figure 8. Cartoon showing the location of the iron absorption 
region. The white dwarf and the disc are not drawn to scale. The 
absorption region is in the outer disc and elevated at a height h 
above the plane of the disc, so as to lie in the line-of-sight between 
the observer and the white dwarf. The intersection of the line-of- 
sight and the gas layer is at a radial distance ri„ from the white 
dwarf. The region of the gas layer which we observe is dependent 
upon the phase of the observation. 

three models in which the gas layer is above the plane of the 
disc. We choose the heights of these three models such that 
ri n is 5, 10 and 15 white dwarf radii from the disc centre. We 
see that when the absorption region is in the plane of the 
disc the model line profile is much broader than is observed. 
If we consider the models in which the absorption region is 
in the outer disc, we see that as the distance of the absorb- 
ing region from the white dwarf increases the profile of the 
line narrows, and loses the broad wings which are observed 
at smaller radii. The models in which ri„ is between 5 and 
10 white dwarf radii are a good match to the observed line, 
suggesting the iron curtain material is at a height h of 1.6 
to 2 .2 white dwa r f radi i above the disc plane. 

iHorne et al.l (1 19941 ) attributed a velocity dispersion of 
60km/s to the iron curtain material they detected in the ul- 
traviolet. In Figure [5] the dashed lines show the same four 
models when we assume the RMS line width to be equal to 
this velocity. The effect is to broaden the line profiles con- 
siderably. Many of the lines seem somewhat narrower than 
all of these models, and so we do not think the velocity 
dispersion of the material we observe is this high. The ve- 
locity dispersion may however be somewhere between this 
value and our initial assumption of lOkm/s, which is a re- 
alistic minimum. This would imply a higher ri n , which in 
turn would place the material closer in height to the plane 
of the disc than we determined for the lOkm/s models. 

3.6 The shape of the accretion disc 

As Figure [7] shows, the radial velocities of the iron curtain 
do not track the velocity of the white dwarf, which we can 
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Figure 9. We show here three representative Fe lines from the 
iron curtain, plotting the averaged line in velocity space, with 
the orbital velocity variations removed. Offset below these, we 
show theoretical line profiles generated using the simple gas layer 
model described in Section I3.5I The models are offset for clarity. 
We expr ess all distances in te rms of the white dwarf radius R\ 
(0.01-Rq . iLittlefair et al.l2008t) . We plot three models in which the 
absorbing layer is above the plane of the white dwarf. We choose 
the height in these three layers so that the minimum radius of 
the observed region of the gas layer ri n is at a distance of 5, 10 
and 15 white dwarf radii from the white dwarf itself. We also 
plot a model in which the gas layer is in the plane of the white 
dwarf, in which case all of the absorbed light originates in the 
lower hemisphere of the white dwarf, and the absorption occurs 
in the inner disc. Finally, for each model we plot two lines, one 
for an assumed RMS line width of lOkm/s (solid), and on for an 
assumed width of 60km/s (dashed). 



predict with confidence given OY Car's well-determined pa- 
rameters and ephemeris. The difference between the iron 
curtain and white dwarf velocities is evidence of significant 
radial motion of the material causing the absorption. For 
instance at phase 0, we see material receding from us at 
150km/s while the white dwarf is moving at 66km/s away 
from us (the systemic velocity). This outer disc material, 
which is on our line of sight to the white dwarf, thus has a 
radial component of 84km/s towards the white dwarf. Tak- 
ing this reasoning further, if we assume that the material we 
see at all phases is on a single path around the white dwarf, 
we can integrate the radial velocity to find its radius as a 
function of azimuth using 



R = Ro + 



Vr 



where Ro is an arbitrary initial radius. The sini projection 
factor allows for our assumption that the motion we see is 
parallel to the orbital plane. We can convert this into an 
integral over azimuthal angle 9 using uj — d9/dt, whereby 



R=Ro + 



Vr 



d0'. 



We take uj to be given by the standard Keplerian relation for 
a fixed radius Ro. This ignores variations that must occur 
given the varying radius, but should be good to first order. 
We select Rq to match the observed outer disk velocity. The 
azimuthal angle is directly related to the phase <j> of Figure 
[7]by 9 = -2tt0. 

To enable the integration, we fitted the velocities with 
four Fourier components. The result is plotted in the left- 
hand panel of Figure 1101 A small but significant distortion 
of the outer disc is evident with the main asymmetry be- 
ing in the {/-direction. T his stands i n cont rast to three-body 
orbits such as those of IPaczvnskil (1 19771 ) which are asym- 
metrical in x but not y. If the distortion deduced from the 
iron curtain is real, then we might expect to see it reflected 
in the Doppler maps, assuming that the emission lines orig- 
inate from associated gas. We therefore also predicted the 
equivalent path in the Doppler image shown in the right- 
hand side of Figure [TOl (dotted line). Due to the v oc R~ x ^ 2 
relation between orbital speed and radius, this is even less 
distorted than the real-space figure, but we believe that it is 
in rough accord with the image, although the latter appears 
more distorted. We show this by plotting points marking 
the maximum of the emission ring. It can be seen that at 
most azimuths these points lie closer to the dotted, distorted 
path than they do to the solid, symmetrical circular path. 
We emphasize that the agreement here is between two inde- 
pendent lines of evidence: the path of the outer disc deduced 
from measurements of the iron curtain material as it is seen 
in absorption against the white dwarf versus the emission 
ring seen in the Ca II. The disc is shown to be non-circular, 
and the distortion although small, is enough to render any 
measurement of Ki from the emission lines entirely useless, 
as is already strongly suggested by Figure [7] It will be inter- 
esting to establish whether hydrodynamic simulations can 
match our observations. We cannot tell given our single orbit 
whether or not the distortion we see is fixed in the rotating 
frame of the binary or not. Further epochs of observations 
will be needed to determine this. This distortion could, for 
example be the disc ellipticity which is observed as super- 
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Figure 10. Left: The shape of the accretion disc, determined by fitting the iron curtain lines. The x and y coordinates are scaled by the 
binary separation a. When compared to a circle centred on the centre of mass(solid line), there is an distortion of the outer disc (dotted 
line) with the main asymmetry being in the y-direction. The dashed line shows the extent of the Roche lobe. Right: Here we transform 
the positional information from the iron curtain fit into velocity coordinates. Again, the dotted line shows the distortion of the outer 
disc as determined from the iron lines. We plot also the Ca II doppler map, and the large dots denote the maxima of the emission ring. 
These maxima appear to track the iron curtain velocities more closely than the symmetrical circular path (solid). The two points in the 
centre of the plot show the centres of mass of the system and the white dwarf. 



humps in the light curves of SU UMa-type CVs (such as 
OY C ar) in the aftermath of superoutbursts in these sys- 
tems (He ssman et al.lll992T ). A superoutburst was observed 
in OY Car from 2009 November 7 to November 210, three 
months prior to our X-shooter observations. 



4 CONCLUSIONS 

In this paper we have presented high resolution spectra of 
the eclipsing cataclysmic variable OY Car. The primary mo- 
tivation of this work was to obtain a spectroscopic determi- 
nation of the donor star orbital velocity, to compare with 
the value which has been inferred from photometric models 
of the eclipse features (K-z = 470 ± 7km/s). We find two K I 
doublets in the J-band track the motion of the donor star, 
and subsequently determine K2 = 470.0 ± 2.7km/s, which is 
consistent with the photometric determination. As well as 
the absorption features from the donor star and the usual 
emission features seen in cataclysmic variables, we also de- 
tect a forest of absorption lines concentrated in the 3100 - 
3500A region. We identify the the majority of these features 
to be Fe I and Fe II lines, and they have an approximately 
sinusoidal radial velocity curve with a phase and amplitude 
which is not consistent with an white dwarf or donor star 
origin. We attribute these features to the 'iron curtain' of ab- 
sorbing material from the outer disc, which ve ils the white 
dwarf in this system and was first reported bv lHorne et al.l 
1 19941 ). The velocities of these lines do not track the white 
dwarf, reflecting a distortion of the outer disc that we see 
also in the Doppler images. This is evidence for considerable 

1 Determined from data obtained by the American Association 
of Variable Star Observers jhttp:/ /www. aavso.org/j l 



radial motion in the outer disk, at up to 90km/s towards and 
away from the white dwarf. 
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